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Abstract 

Resonant field amplification, which may lead to locking or error field penetration, involves driving 
a weakly stable tearing perturbation in a rotating toroidal plasma. In this paper it is shown that 
the locking characteristics for modes with finite real frequencies are quite different from the 
conventional results. A calculation of the tearing mode amplitude assuming modes with frequencies 
in the plasma frame shows that it is maximized when the frequency of the stable backward 
propagating mode {—uor) is zero, i.e. when v = uor/k. Even more importantly, the locking torque 
is exactly zero at the mode phase velocity, with a pronounced peak at just higher rotation, leading 
to a locked state with plasma velocity v just above the mode phase velocity in the lab frame. Real 
frequencies are known to occur due to the Glasser effect [H [2] for modes in the resistive-inertial 
(RI) regime. This therefore leads to locking of the plasma velocity to just above the phase velocity. 
It is also shown here that real frequencies occur over a wide range of parameters in the visco- 
resistive (VR) regime with pressure, and the locking torque is similar to that in the RI regime. 
The same is true in other tearing regimes, with real frequencies due to diamagnetic propagation. 
Other nonlinear effects and the possibility of applying external fields of different helicities to drive 
sheared flows in toroidal plasmas are discussed. 
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I. INTRODUCTION 


Resonant field amplification [7] (RFA) is the response of a magnetically confined toroidal 
plasma to an externally applied nonaxisymmetric magnetic field, an error field. Error fields 
are generally unavoidable, and in some cases such fields can be applied purposefully, as 
discussed below. RFA can lead to error field penetration (EFP) (EFP)|3]-[6l, [20l [ 2 T] or 
locking. The accepted view is that EFP involves locking the plasma rotation to just above 
zero velocity in response to the error field. If this occurs it can cause a loss of the benefit of 
sheared rotation on plasma stability, and on occasions can lead to disruptions in tokamaks. 
It is well known that RFA occurs by means of coupling to a stable tearing mode, and that 
the response is maximized when the tearing mode is near marginal stability [3H5l [71 [ 20 l| 2 T]. 
The dispersion relation for a spontaneous mode (with 'ip{ryj) = 0 ) is A( 7 d) = A', where 
7 d is the growth rate 7 + ikv in the plasma frame and A' and A are the usual constant-V' 
matching parameters from the tearing mode outer and inner regions, respectively. Here, v 
is the toroidal plasma velocity (E x B plus parallel velocity) at the mode rational surface. 
The steady-state magnetic perturbation at the edge of the tearing layer at the mode rational 
surface {k\\ = 0 ) is called the reconnected flux in the context of the constant-?/^ approximation, 
and is known to have the form 


= -7U- 


l 


21 




( 1 ) 


A' — A(ikv) 

Here, 'fi{ryj) is the amplitude of the error field applied at the wall, is the radius of the mode 
rational surface, and I 21 is an inductance coefficient. We describe the results here in terms 
of large aspect ratio cylindrical geometry, although one toroidal effect, namefy favorable 
average curvature for safety factor g > 1, is of importance and is included. The quantities 
A, A' and /21 are normalized so that the plasma radius is equal to unity. 

For tearing modes in the viscoresistive (VR) regime [H |9], with resistivity r} and ion 
viscosity fi, A ( 7 ) has the form 


= ^5/6 1^/^11/351/3^*^ = ( 2 ) 

where is a dimensionless positive constant of order unity. From Eq. 0 it is seen that 
spontaneous VR tearing modes have real 7 for all A' and Eq. 0 shows that the largest 
response to 'fi{ryj) for stable plasmas is near marginal stability, where \A.' — ikvTyr \ is smallest. 
This occurs for A' ^ 0— with —u)r + fcn = 0, i.e. n = 0. 
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The toroidal torque on the plasma, which comes from the Maxwell stress in the non-ideal 
layer, takes the quasilinear form [HI 0] 

iVm =(3) 

(This quantity is the integral over the tearing layer of j\\Br and is formally a force in this 
reduced model. In a tokamak, because of strong damping of the poloidal flow, the toroidal 
force component is most important, which gets multiplied by the factor rt/q{rt)R, where q is 
the safety factor and R is the toroidal major radius. Thus there is a direct relation between 
Nm and torque along the tokamak symmetry axis from toroidal force.) Using Eq. ([^ it is 
readily seen that this torque in the VR regime is 

^ ... 

- 2 A'2 + fcVT2/ ^ ’ 

behaving as —v/{y\ + with extrema at k‘^v^ = A'^/r^^, and symmetry about v = 0, 
Nm{—v) = —Nm{v). The quantities \'ip{rt)\‘^ and as a function of h = kvTyr are shown 
in Fig. illustrating the behavior around u = 0. 

Based on these results, it has been shown [HHSl C] that models with a viscous response 
torque at the layer Ny = — v) (where Nq represents a balance between the momentum 

source and viscous drag, with no error held to give an equilibrium flow Vq) can have three 
rotation states, intersections with N.m-\-Nd = 0, for intermediate values of |'^(ru,)|, as shown 
in Fig. The bifurcation diagram illustrating these states is shown in Fig.[^ The fastest is 
the rapidly rotating or high slip state, and the slowest is the penetrated or locked state. The 
middle intersection represents an unstable or forbidden rotation. For small |^/;(r^)| or large 
Nq (or large Vq), only the rapidly rotating state exists; for increasing |^/;(r^)| or decreasing 
Nq there is a bifurcation at 'tpyji past which the locked and intermediate states can exist. 
For further increasing or decreasing Nq, there is a bifurcation at 4’w2 past which 

the two upper states coalesce and only the locked state can exist. Since '^u )2 > these 
bifurcations have hysteresis [HI IH HDl (HI 120] • 

II. RESISTIVE-INERTIAL (RI) REGIME 

In the resistive-inertial (RI) regime, the quantity A takes the form 

^ ' = 7 ^ U/ , (5) 


3 






- 3 - 2-10 1 2 3 

V 


FIG. 1. Plot (a) of reconnected flux magnitude vs. v = kvTyr for a weakly stable VR mode. 

In (b) are shown the torque —Nm vs. v as well as the viscous torque curve Ny ex vq — v. In this 
parameter range there are three equilibria, with the locked equilibrium having very small f) ^ 0. 





FIG. 2. Bifurcation diagram for the VR regime, a plot of the the equilibrium flow vs. the wall 
perturbation amplitude. 
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where is another positive dimensionless constant. In the presence of equilibrium pressure 
gradient p'{rt) and field fine curvature at the tearing layer, parallel dynamics, perpendicular 
compression and particle transport, it is known that for sufiiciently high plasma P this is 
modified to give|8l [l2] 

A, ^ (A - S/73/2), ( 6 ) 

where A and B oc Dg are constant and Dg = — 2 rp'/BgR'^q'‘^ is the Suydam parameter. For 
Dg > 0(S > 0), this result leads to the observation that tearing modes are destabilized for all 
values of A' [HlIE], becoming electrostatic with 7 r\j (Dg/A'p as A' ^ —00 [ 12 ]. In toroidal 
geometry Dg is replaced by the Mercier parameter D = — 2 rp'(l — (p)/BlB?q'‘^ = — 

giving S < 0. Again, q is the safety factor and R is the toroidal major radius. The 1 — q^ 
factor is due to the average curvature in toroidal geometry and leads to D < 0 for q{rt) > 1 
and p'{rt) < 0. 

For favorable curvature D < 0 the dispersion relation A(7) = A' gives two real roots 
with 7 > 0 for A' > A^m oc \DpA_ Por A' < Amin, there are complex conjugate roots 
with 7 = 7r T in fhe plasma frame; these roots have 7r < 0 (are stable) if A' < Ac; the 
quantity Ac is also proportional to |Zi>| 3 /®.[Il [2| 


The possibility of complex conjugate roots (roots for the spontaneous modes with real 
frequency) near marginal stability 7^ < 0 or A' < Ac, leads to the main result in this paper. 
The reconnected flux from Eq. 0 satisfies 
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hi 


\A'— A(ikv)P 

and combining this and Eq. Q we find 

2 {A'— Ar{ikv)y + Apikvy 


(7) 


The denominators |A' — A{ikv)p are minimized {\ip{rt)P is maximized) when the roots 7 
corresponding to the spontaneous mode, with A' = A{'jci), are closest to zero in the complex 
plane. This occurs when Aojr + kv is closest to zero, i.e. the frequency of the backward 
wave in the lab frame —Ur + kv is zero. See Eig. This allows maximum interaction with 
the zero frequency error field. Note that at this point Apikv) is approximately zero so, by 
Eq. Q, Nm is very close to zero where v = ujr/k. 
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FIG. 3 . Sketch of the locus of roots in the RI regime with favorable curvature. The blue arrows 
indicate the change in the roots as v increases, with the backward wave complex frequency coming 
closest to zero for v = uir/k. 

The quantities and are shown in Fig. the latter with the viscous torque 

Ny{v). Note that the reconnected flux is peaked away from u = 0 , approximately where 
V = ujr/k. Note also that increases rapidly for v > ujr/k, so that the error fields tend 
to lock the plasma to just above the phase velocity of the tearing mode uJr/k, as shown in 
Fig-i rather than to just above zero plasma velocity as in Fig. This effect produces a 
positive asymptotic value of plasma rotation v uJr/k for the locked state as boundary 
held is increased, as shown in Fig. The upper line shows three equilibria, but with the 
leftmost intersection having v > uJr/k rather than u > 0 . For spontaneous modes which 
are very weakly damped and have appreciable phase velocities, the locked state can have 
significant plasma rotation. Also note that the intersection with the lower line also indicates 
three equilibria, but with two negative values of u, and one with v > uJr/k. All of these 
locked states with finite velocity exist in spite of the fact that the magnetic field driven by 
the error held is steady, i.e. locked to the error held. This effect is related to the fact that 
for RI (but not VR) resistive wall tearing modes, small rotation can destabilize the modes 
by maximizing the coupling of the tearing modes to the resistive wall| 13 |. even for Dg = 0 . 

We have formulated the inner layer equations in the RI regime with p', curvature and par¬ 
allel dynamics, but ignoring perpendicular compression and perpendicular particle transport, 
in order to determine the importance of the last two effects on the stability and propagation 
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FIG. 4. Plot (a) of |'0(rt)p vs. v = kvTri for a weakly stable RI mode with the Glasser effect, 
showing peaks at 'u = ujr/k^ i.e. v — uOrTri- In (b) are shown the torque curve —N^ vs. i), with 
Nr^ = 0 close to the phase velocity v — uJrTri = 0.12; also shown are two lines Ny = iVo('^o ~ '^) 
with different values of N^. Note the differences between this figure and Fig|T] 



FIG. 5. Bifurcation diagram for the RI regime, indicating the locked (penetrated) state, the 
rapidly rotating state (high slip state) and the middle unstable state (labelled ’forbidden’.) Note 
the positive asymptote v Ur/k as increases. 
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of the spontaneous modes. Using a scaling specific to this constant-^ regime, we obtain 


d^W 


-i^w 


GQ^ 


g2 + 62 ^: 


:1U = - 1 + 


Gb‘^ 


g2 + 62^2 




( 8 ) 


where 0 = -iaS^xpoW(^)/pot], r-n = S^,S = G = - 2 m?Bl{rt)p'/Bl 

Dg, and Q = ^r^i. Also, we have G = Gq/Q^G and 6^ = Xhe effect of toroidicity 

is included by taking Go —> (1 — h(u)^)G^o as above, leading to Go < 0 for p'{rt) < 0 and 
q{rt) > 1 . The constant-'^ small parameter is e = hA' ~ yv'e assume p' ~ e, leading to 
G ^ b ^ 1. The matching to the outer region gives 


A' = 7 r (1 _ ^U') d? = ^A.(Q) = 7=/‘‘tXa.(Q) = 0='‘‘A,(Q) = A(Q). 

d J-oo V 


(9) 


The factors l/(g^ + b^^^) oc 1/(7^ + in Eq. (jS) represent sound wave propagation 

parallel to B; the factor Q"^/(Q^ + 6^^^) on the left is stabilizing; the term oc G on the right 
is destabilizing but diminished by the l/(g^ + factor. 

Equation (j^ and therefore As(g) are invariant under the symmetry Go ^ AGo, bo — >■ 
A^/^6o; Q —t This symmetry has two invariant quantities Gq/^q 0^ p'/Tp and 

G = GoIQ^G, Eor Go positive (unfavorable curvature) and sufficiently large, the As(g) 
curve on the real Q axis has poles (not shown) corresponding to localized unstable electro¬ 
static resistive interchanges Qn. These modes are stabilized {Qn 0 ) for sufficiently large 
sound speed (Go/60 small enough.) Past this threshold we observe numerically the relation 
As(g) ~ A — B/Q^G^ which shows 


As{Q) = Ag{Q = 00 ) 


Go f Go 

WG 


for some function K. A fit to the numerical data for large sound speed Go/b^ <C 1 gives 

2.77Go 


A3(g) = 2.12- 


g 3/2 


1 


0 . 65 -^ 

^0 


( 10 ) 


This with Eq. § shows the form of Refs. [U E], namely A(g) = AQ^G _ i/ 4 ^ where 

B (X Go oc Dg. The unfavorable curvature (Go > 0 ) results of Ref. [12] are recovered. In 
particular, these results show that Q r\j (Go/|A'|)^ as A' — 00 . 


The fit in Eq. (10) works well in the favorable curvature case (Go < 0) also, again for 
sufficiently large sound speed (small |Go/6o|.) Unlike in the unfavorable curvature case, 
there is no evidence of unstable electrostatic resistive interchanges. That is, these modes 
have complex Q and do not show up as poles in A(g) on or near the real axis. These results 









for favorable curvature show that the behavior observed in Refs. in El holds qualitatively 
using only parallel dynamics, without including the divergence of the E x B drift and 
perpendicular particle transport. In particular, there are complex roots for A' < Amin = 
min(A((5)) ~ and these roots become stable for A' < Acrit ~ The locus of 

roots is shown in Fig. 

The quantities , —Nm and as functions of the plasma rotation rate v shown in 

Fig. 1^ for A' < Ac can be compared with those in Fig. [1} For such values of A' the mode 
is weakly stable and the torque is peaked just to the right of v = ojr/k, and negative to the 
left of this rotation value. The slope of Nm{v) is steep near v = u)r/k, so the locked state is 
just to the right of the phase velocity. As in the VR case in Fig. (with Ur = 0 ), there can 
be two other roots, the rapidly rotating state and an unstable middle root. 

It is instructive to consider the case with p'{rt) = 0 (Go = 0 ) in the RI regime. As 
shown in Fig. the reconnected flux |'?A(rf)p in the RI regime with and A' < 0 also has two 
peaks at kvTri 7 ^ 0, corresponding to stable complex conjugate roots with Q r\j |A'|4/5e±4W5_ 
However, the torque curve is not qualitatively different from the VR curve in Fig. for zero 
pressure gradient, since Nm ~ —Im(iA:n)^/4 — —sgn(/cn)|A:np/4 sin( 57 r/ 8 ), which goes to zero 
only at n = 0. It can easily be seen that the simple form for the VR regime (without pressure 
gradient) Nm ~ —v/{vf + n^)( as in Refs. P-OC]) does not hold quantitatively. 


III. VISCORESISTIVE (VR) REGIME 


For the VR regime, with pressure gradient, curvature and parallel dynamics, but without 
perpendicular compression or particle transport, the inner region equation for the stream- 
function W takes the form 


d^W 


^^IF- 


GQ2 


W = 


^^4 ^ Q^ + b^e 

where r — Vt = with S = Q = 

b = aScgTyr/B. The matching condition is 

A' = ^ r (1 - ^IF) = ^A^Q) = 

V J-00 V 


G (X Dg is defined as in the RI case, and 

^TnrAsiQ) = QAsiQ) = A(Q). (12) 


We also have 


Q' 
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FIG. 6. The reconnected flux magnitude (a) \'ip{rt)\‘^ in the RI regime with Go = 0 shows peaks 
at flnite v = kvTri, corresponding to complex roots, but the torque shown in (b) goes to zero 
only at u = 0, so locking to a nonzero velocity v does not occur. 


and 

b = bo. 

Equation 0 has the symmetry Go ^ AGq, bo —>■ Xbo, Q XQ, so the invariants are 
Go/bo and Gq/Q. Notice that this symmetry extends to A < 0, i.e. Q —>■ —Q, Go —Go 
(and trivially to bo —bo-) The quantity is unchanged, but A(( 5 ) —t — A(( 5 ). Thus the 
results for unfavorable curvature (Go > 0) can be applied directly to the favorable curvature 
(Go < 0) case. 

We will describe the results first in terms of favorable curvature (Go < 0). For large 
values of \Go/bo\ ~ \p'\l\/p the plot of As(( 5 ) on the positive real Q = v = axis shows 
poles, again corresponding to electrostatic resistive interchanges. As \Go/bo\ decreases most 
of these modes become stabilized, but the last two coalesce at |Go/6o| = 9.6 for positive Q, 
leaving a nonmonotonic A(( 5 ) curve as shown in Fig. This curve is nonmonotonic for the 
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FIG. 7 . Plot of A{Q) in the VR regime with favorable curvature for Go/bo = — 1 . 6 . 



FIG. 8. Locus of roots showing that there is a Glasser effect in the VR regime. 

wide range 1.2 < IGq/^oI < 9 . 6 . For A' just below the minimum on this curve, or just above 
the maximum, there are complex roots. A plot of the locus of these roots is shown in Fig. 
We obtain a reconnected flux curve and a torque curve shown in Fig. similar curves to 
those of in Fig. 1 ^ for the RI regime. (The example in the VR regime in Fig. is for a much 
more weakly damped mode, leading to more peaked functions.) 

For unfavorable curvature (Go > 0), the symmetry Go —Go, Q —Q shows that 
similar complex roots occur, but they are weakly growing when the roots for Go > 0 are 
weakly damped. 


11 









FIG. 9. Plot (a) of reconnected flnx magnitude |'0(rt)p vs. v = kvTyr for a weakly stable VR mode 
with pressure gradient, curvature, parallel dynamics, in the regime with a non-monotonic form of 
A((5). Again, the reconnected flux magnitude is maximum near v = oor/k. In (b) are shown the 
torque curve —N^ vs. v, with = 0 at the phase velocity v = oJrTyr = 0.13 with the two lines 
Ny = Ao(^o “ v) having different values of Nq. Here there can be as many as five equilibria, and 
the intersection at the lower positive v, the locked state, again has v > rather than h > 0. 


IV. CONCLUSIONS 

Our major new result relates to the torque Nm on the plasma as a fnnction of velocity 
V for tearing mode regimes with real frequencies, snch as the resistive-inertial (RI) regime 
with the Glasser effect. We find that N^iy), which is concentrated in the tearing layer of 
the driven mode, is exactly zero where v = ujr/k. This torque Nm{v) is peaked just to 
right oi V = ojr/k, as shown in Figs. and If we model the viscous torque driving the 
plasma rotation as Ny = Nq{vo — v), then we find that an error field can cause locking to 
a plasma velocity v with v > Ur/k. Thns, the magnetic pertnrbation locks to the error 
field at rest in the laboratory frame, bnt the plasma locks to a velocity just above the finite 
phase velocity v > uJr/k of the spontaneous tearing mode. Fnrther, onr results show that 
the plasma velocity asymptotes to this phase velocity as oo or Vq ^ 0 . 
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We have also shown that in the VR regime with pressure gradient and parallel dynamics 
in the tearing layers, there is a realistic and wide range of parameters for which weakly 
stable tearing modes with complex oj ~ zbcj^ + W can occur. That is, there is a Glasser 
effect in the VR regime as well as in the RI regime [U [2]. In both regimes, the presence 
of real frequencies ojr in the plasma frame implies that the driven tearing mode has a very 
peaked maximum amplitude where the plasma velocity v equals the phase velocity ujr/k. 

Tearing modes in the presence of diamagnetic effects can also have real frequencies Ur oc 
(u*, with a reduction in growth rate 7 duni. It is expected that diamagnetic propagation 
will cause a similar locking of a plasma to a rotating state, v > u)r/k. For such modes, 
however, there will be no n —v symmetry because the modes do not occur in complex 
conjugate pairs. 

For sufficiently large error fields, the mode amplitude in the locked state can be high 
enough to be in the Rutherford regime [TH]. where the magnetic island width w is compa¬ 
rable to the tearing layer thickness. In such cases, the linear theory for the mode (and 
quasilinear theory for the torque) should still be qualitatively accurate. For reasonable 
DIII-D parameters with moderate plasma rotation, we find that the locking threshold is 
sufficiently small that the tearing mode amplitude is below or comparable to that for the 
Rutherford regime, so that the above conclusions are valid. 

At even higher error field amplitudes, sound wave propagation can flatten the pressure 
around the island [T^ and this weakens the mode propagation. For DIII-D parameters with 
fast rotation, the locking threshold is for a large error field, so that typically the plasma will 
not lock. However, in the presence of such a large error field, locking by a necessarily large 
error field can lead to the formation of a quite large island, and the finite frequency effects 
discussed in this paper, due to pressure gradient, may not be present. 

We have one caveat about the linear tearing mode regimes and the nonlinear regimes 
discussed above. These are useful for guidance but, for example, when analyzing modes in 
the VR regime, the mode may move into another regime as the velocity increases. Also, the 
correct regime may be different for the locked and unlocked states in the same bifurcation 
diagram. 

These results provide a new viewpoint on the locking phenomenon: the locking force by 
itself leads to a plasma velocity equal to the phase velocity of the spontaneous mode. (The 
actual velocity of the plasma must be calculated using a balance the other torques, a mo- 
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mentum source due to neutral beams.) As an example of an application of this framework, 
our results suggest the possibility of applying error fields with a spread in Fourier spectrum. 
Then locking may occur in each tearing layer to just above the phase velocity for its Fourier 
harmonic. If these phase velocities are not equal, then viscosity (e.g. NTV) should lead to a 
sheared rotation between these mode rational surfaces. This possibility might be useful for 
maintaining fiow shear to help stabilize other modes. 
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